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ABSTRACT 

We present an analysis of the effects of beam deconvolution on noise properties in CMB measurements. The analysis is built around 
the artDeco beam deconvolver code. We derive a low-resolution noise covariance matrix that describes the residual noise in decon¬ 
volution products, both in harmonic and pixel space. The matrix models the residual correlated noise that remains in time-ordered 
data after destriping, and the effect of deconvolution on this noise. To validate the results, we generate noise simulations that mimic 
the data from the Planck LFI instrument. A^ 2 test for the full 70 GHz covariance in multipole range £ = 0-50 yields a mean reduced 
X 2 of 1.0037. We compare two destriping options, full and independent destriping, when deconvolving subsets of available data. Full 
destriping leaves substantially less residual noise, but leaves data sets intercorrelated. We also derive a white noise covariance matrix 
that provides an approximation of the full noise at high multipoles, and study the properties on high-resolution noise in pixel space 
through simulations. 
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1. Introduction 

Several methods have been proposed for deconvolution map¬ 
making for CMB (cosmic microwave background) measure¬ 
ments (|Arrmtage & Wandelt|2004| |Armitage-Caplan & Wandelt 
|2009[ [Harrison et al.||201 1 [ |Keihknen & Reinecke||2012| l. The 
aim of deconvolution mapmaking is to produce a map that is 
free from effects of beam asymmetry. For instance, the image 
of a point source in a deconvolved map appears symmetric, in¬ 
stead of being elongated according to the beam shape. More im¬ 
portantly from the point of view of cosmology, deconvolution 
mapmaking produces a map with a unique effective beam win¬ 
dow that takes the same form in all map pixels. This happens at 
the cost of more complicated noise structure. A simple binning 
operation keeps white noise uncorrelated: white noise in the in¬ 
put time-ordered information (TOI) translates into uncorrelated 
noise in the binned map. This is not the case for deconvolved 
maps. The deconvolution operation creates correlations between 
pixels. Correlated noise present already in the input TOI compli¬ 
cates the noise structure further. 

In this work we study the impact of beam deconvolution on 
noise properties in CMB maps. We build our treatment around 
one particular deconvolution code, artDeco. The artDeco code 
( Keihanen & Reinecke|2012) was developed for beam deconvo¬ 
lution of absolute CMB measurements. The code takes as in¬ 
put the time-ordered data stream along with pointing informa¬ 
tion and a beam description, and produces as primary output 
the harmonic coefficients axt m (X = T, E, B ) of the sky. From 
these coefficients one can further construct a sky map that is free 
from beam asymmetry effects. The computation is carried out in 
harmonic space. In doing so the method breaks the assumption 


shared by pixel-based mapmaking methods that the sky signal is 
constant within a pixel. 

We investigate two kinds of input noise: pure white noise 
and destriped 1 // noise. The starting point in the latter is a sit¬ 
uation where TOI is first destriped with a general destriper with 
noise prior. One source of such data is the Madam map-making 
code ( |Keihanen et al.||2005| |2010| ), which we use in this work. 
The code subtracts an estimate of the correlated noise compo¬ 
nent from the TOI and yields a cleaned TOI stream where resid¬ 
ual noise is dominated by white noise. This cleaned TOI is then 
provided as input to the artDeco beam deconvolver. 

The noise in the output map is correlated for two reasons: 

- There is residual correlated noise in the TOI stream that can¬ 
not be entirely removed by destriping techniques. 

- The deconvolution process itself changes the noise proper¬ 
ties. 

We aim at deriving a noise covariance matrix (NCVM) that cap¬ 
tures both phenomena. 

We validate our results using simulated Planck LFI data. 
However, we aim at keeping the method more general. The re¬ 
quired conditions are as follows: 

- The experiment records the absolute signal, not differential. 

- The experiment offers full or nearly full sky coverage. 

- Beam shapes are known. 

- Sky sampling is dense compared to beam size. 

- Noise is piecewise stationary and its spectrum is known. 

- The data is cleaned of correlated noise through destriping or 
another method with equivalent results. 

We do not require a particular scanning pattern, or a particular 
noise spectrum. 
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Our calculations are carried out in harmonic space, which is 
the natural domain for beam deconvolution. We construct a har¬ 
monic noise covariance matrix that describes the noise proper¬ 
ties of the harmonic axtm coefficients, which are the primary out¬ 
put of artDeco. Because the problem is huge, it is not possible 
to construct the covariance matrix for all multipoles. In practice, 
the noise covariance matrix approach is limited to multipoles 
i % 100. 


We further produce a pixel noise covariance matrix, which 
describes the noise properties of deconvolved sky maps, con¬ 
structed from the cixi m coefficients through spherical harmonic 
transform. A noise covariance matrix is an essential ingredient in 
many CMB power spectrum estimation methods. A method for 
producing a noise covariance matrix for un-deconvolved maps 
has been presented in Keskitalo et al. (20101. Its application to 
Planck LFI data is presented in Planck Collaboration VI ( [2015) 1. 

The paper is structured as follows. In Sect. [2] we perform a 
series of full-resolution simulations, to gain a general view of 
the impact of deconvolution on residual noise. In Sect. [3] we fo¬ 
cus on white noise. We derive a white noise covariance matrix, 
and compare its predictions with Monte Carlo simulations. The 
main result of this paper is the derivation of a full low-resolution 
noise covariance matrix, presented in Sect. [4] In Sect.[5]we per¬ 
form a thorough validation of the harmonic covariance matrix 
through Monte Carlo simulations. We show results from x 2 tests 
and from noise bias comparison. In Sect. [6] we perform analysis 
in pixel space. We summarise our results in Sect. [ 7 ] 


2. Simulations 

We start by performing a series of full-resolution noise simu¬ 
lations to assess the impact of deconvolution. Our simulations 
mimic the data from the LFI instrument of Planck experiment. 
We created a time-ordered data stream, which we filled with sim¬ 
ulated noise, and fed it as input to the artDeco deconvolver. The 
simulations were run on Sisu, the Cray XC40 supercomputer of 
CSC, Finland. 


2.1. Generating noise 


We regenerated the LFI detector pointing with Levels software 
(Reinecke et al .|2006 1. The simulated mission covers 4 years of 
data. We used the internal noise generator of the Madam map¬ 
making code to generate a noise time stream. The noise proper¬ 
ties mimic those of the real Planck data. We considered the three 
LFI channels: 30 GHz, 44 GHz, and 70 GHz. The noise param¬ 
eters used in simulations were taken from iPlanck Collaboration! 
II ( 2015| >. For convenience, the parameters are listed in Table[I| 

We consider two types of noise. In the first case we generated 
pure white noise, with variances taken from Table |T] The vari¬ 
ance is assumed constant in time, but varies between radiome¬ 
ters. In the second series of simulations we generated realistic 
1 // noise, again with parameters from Table[l] We destriped the 
noise stream with Madam, and stored the cleaned data on disk. 
The data was then fed as input to artDeco. 

The destriping procedure follows the one applied to LFI 
mapmaking (Planck Collaboration Vl]|2015| >. We applied the 
same destriping mask, and used the same HEALPix map resolu¬ 
tion A sid e=1024 (3.44’), and same baseline length (Is for 70 GHz 
and 0.25 s for 30 GHz) as in real Planck mapmaking. We refer to 
these simulations as realistic noise simulations, in contrast to the 
white noise simulations. In both simulations we applied realis¬ 
tic Planck flags to the TOI stream. This excludes roughly 8% of 


Table 1 . Noise parameters used in simulations: Knee fre¬ 


quency (/knee), slope (/3), and white noise rms (cr). From Planck 
[Collaboration II| ( |2015| l. 


/alee / 


Horn 

M 

S 

M 

S 

M 

S 

70 GHz 

LFI-18 

14.82 

17.78 

-1.060 

-1.180 

4.553 

4.146 

LFI-19 

11.72 

13.71 

-1.207 

- 1.111 

5.144 

4.926 

LFI-20 

7.96 

5.67 

-1.198 

-1.298 

5.212 

5.507 

LFI-21 

37.89 

13.27 

-1.247 

-1.205 

4.003 

4.971 

LFI-22 

9.68 

14.80 

-1.424 

-1.237 

4.356 

4.715 

LFI-23 

29.65 

59.03 

-1.073 

-1.211 

4.476 

4.790 

44 GHz 

LFI-24 

26.78 

88.30 

-0.942 

-0.908 

3.159 

2.734 

LFI-25 

20.07 

46.37 

-0.845 

-0.904 

2.834 

2.698 

LFI-26 

64.42 

68.19 

-0.918 

-0.758 

3.295 

2.887 

30 GHz 

LFI-27 

174.53 

108.79 

-0.927 

-0.907 

1.605 

1.729 

LFI-28 

130.08 

43.08 

-0.931 

-0.900 

1.812 

1.633 


the data. Madam has the capacity of storing the destriped TOI on 
disk in “4D map” format. The four dimensions in this context re¬ 
fer to the three angles 8, <p, it/ which define the detector pointing, 
and pointing period index as fourth dimension. Pointing period 
refers to a period of typical duration of 40 min, where the Planck 
scanning pattern follows a fixed ring on the sky. 

Madam generates a three-dimensional grid by discretising the 
detector pointing angles. The pixelization obeys HEALPix pix- 
elization in 8 and 4>, and the 4/ angle is divided uniformly into 
Vp si bins. Another parameter ;V SK | C controls the resolution in 8, <p. 
The TOI is binned on the grid according to its pointing, and 
the grid is stored by pointing period. The data volume is com¬ 
pressed typically by a factor of 20 compared to the full TOI. 
Radiometer-specific parameters required in subsequent analysis 
steps are stored as keywords in the 4D object header. For the 
present study we used resolution A/id e =1024, Ap S i=4096. At this 
resolution the 4D maps take roughly 8 GB of disk space per de¬ 
tector. 

In white noise simulations, we use Madam only to generate 
the noise, to apply the flags, and to store the data stream in 4D 
map format. No destriping is applied. 

The 4D map is an intermediate data object that is useful for 
many purposes. In particular, it can further be compressed into 
the “3D map” format that serves as input to artDeco. This step 
is trivial, it consists of coadding all pointing periods on a sin¬ 
gle grid. Further, we can construct a binned map from the same 
4D map input. The binning operation represents “normal” map¬ 
making without beam deconvolution. The combined operation 
of destriping, 4D map production, and map binning, is nearly 
equivalent to the operation of producing a destriped map directly 
by Madam. The only difference is the discretization of the if/ an¬ 
gle, which at the chosen resolution V ps j=4096 has a negligible 
effect on the final map. 

To reduce the scatter in results, we produced several noise 
realizations for each case studied. At 30 GHz, we generated 40 
realizations for each case. At 44 GHz and 70 GHz, where the 
simulations are more demanding, we produced 10 realizations. 
The spectra shown are averaged over available realizations. A 
complete list of simulations is given in Table [ 5 ] 
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2.2. Destriping options 

The full simulated TOI consists of four years of data for all LFI 
radiometers. In cases where we deconvolve only a subset of data 
from a frequency channel, we have two options for the destriping 
step. Either we can use the full data set for destriping, and then 
pass a subset of the cleaned TOI to deconvolution, or we can use 
the same data subset in both processing steps. We refer to the 
two options as full destriping and independent destriping, re¬ 
spectively. Full destriping leads to a smaller residual noise level 
in the final data products, as there is more information available 
for the destriping solution and removal of correlated noise. On 
the other hand, independent destriping has the benefit over full 
destriping that the residual noise in the final products is uncorre¬ 
lated from one data set to another. 

We compared the two destriping options for horn pair 18/23. 
This represents one third of the full 70 GHz data set. The horn 
pair data set was selected for closer examination on the basis 
of low-resolution simulations (Sect. [5j, where it showed as the 
“worst case” subset. 

2.3. Binned maps 

We compare the deconvolved spectrum to the binned map spec¬ 
trum, which we constructed as follows. We bin the same de- 
striped TOI into a sky map as 

m = (P T C; 1 P)- 1 PC n - 1 y (1) 

Here P is the pointing matrix, y is the input time-ordered data 
(TOI), and C„ is a diagonal white noise covariance. If several 
detectors are involved, their data are appended to a single TOI 
stream. 

Binning represents the usual mapmaking process, which 
does not attempt to correct for the beam shapes. We then com¬ 
pute the pseudo spectrum of the map with the Anafast tool of 
the HEALPix package. We refer to the spectrum obtained this 
way as map spectrum. A comparison of the map spectrum and 
deconvolved spectrum gives insight to the impact of the decon¬ 
volution process on noise. We perform the map-binning opera¬ 
tion starting from the same 4D map objects that serve as input to 
the deconvolution procedure. That way we can be sure that we 
have exactly the same inputs in both procedures. 

We deviate from the mapmaking procedure applied in actual 
Planck analysis in one aspect. LFI mapmaking applies the horn- 
uniform radiometer weighting scheme, when combining several 
radiometers into one map (Planck Collaboration VI 2015[ l. Horn- 
uniform weighting reduces leakage of temperature signal to po¬ 
larization through beam shape mismatch, at the cost of slightly 
increased noise compared to the other option, noise weighting. 
In noise weighting, the radiometer weights are constructed as 
wj = l/cr 2 j where crj is the white noise rms for radiometer j, 
from Table |T| In horn-uniform weighting the weights are made 
equal for the M and S radiometers of same horn. Because the 
beam shapes are explicitly accounted for in deconvolution, horn- 
uniform weighting provides no benefit over noise weighting, and 
artDeco internally applies noise weighting. Since our aim is to 
study the effects of beam deconvolution, and not those of dif¬ 
ferent radiometer weighting schemes, we build also the binned 
maps with noise weighting. Horn-uniform weighting is, how¬ 
ever, applied in the destriping phase. This way we have two map¬ 
making methods that take as input the same exact destriped TOI, 
and differ only in the way they deal (or do not) with the beam. 
The detector weighting scheme has a small effect on the residual 
white noise level. We return to this in Sect. 13.31 


2.4. Deconvolution 


ArtDeco yields as output the harmonic coefficients ajc m - UEtm, 
a-Btm- with l = 0...fjnax- Here f max is an input parameter that 
defines the multipole range considered. Another input parameter 
k max sets the level of beam asymmetry taken into account. The 
harmonic beam expansion is constructed for b s a, where 0 < C < 
f max and \k\ < k max . The maximum ( max that can be reached is 
dependent on the beam width. For the full-resolution simulations 
we chose 7 nmx =700 for 30 GHz, 4^=1000 for 44 GHz, and 
(rnax = 1500 for 70 GHz. In all cases we used k max =6. 

The artDeco deconvolution operation can formally be writ¬ 
ten as 

a = (A'C 1 / 1 A) _1 A t C n 1 y (2) 

The solution is built on the assumption that the noise in the input 
TOI is white, or at least dominated by white noise. Matrix A is 
given by (Keihanen & Reinecke |2012j ), 

A L = T,K a D e : k (ojj), o) 

k 


where b is the harmonic representation of the beam, D is a 
Wigner function, which takes as argument the pointing angles 
6, 4>, and \p, and coj is a short for the triplet of pointing angles for 
sample j. 

The spin index s takes values 0 (temperature) and +2 (po¬ 
larization). Index k takes values from 0 to ±k max . Indices t, m 
are the the usual indices of spherical harmonics, £ being related 
to angular scale and m to orientation. In this notation, the out¬ 
put harmonic coefficients are written as a s ( m and the input TOI 
stream as yj. The E, B coefficients are further constructed as 

UElm = ~2 ( a 2{m+a^2lm) ( 4 ) 

UBCm = -g-.(U2fm ~ «-2fm) 

2 1 

In this work we are interested in the impact of deconvolu¬ 
tion on noise. Since both destriping and deconvolution are linear 
processes, we can examine the signal and noise components in 
isolation. From now on we assume that the input TOI consists of 
pure noise (white and correlated). We compute the deconvolved 
spectrum from the output axe m coefficients as 


CxYC - 


1 


2^+1 


z 

m——t 


a Xtm a Yfm 


( 5 ) 


where X, Y stand for T, E, B. We further average the spectrum 
over the available noise realizations. 


2.5. Beams 


Throughout this paper we use Planck LFI beams ( jPlanckl 
Collaboration IV pOB ). We included in the deconvolution pro¬ 
cess the main and intermediate beam components. The contri¬ 
bution of far sidelobes is assumed to be removed before the de¬ 
striping step (Planck Collaboration II 2015 i. 

Table [2] lists some characteristics of the beams per fre¬ 
quency. The FWHM and ellipticity values are based on effective 
FEBeCoP beams (|Mitra et al.|2011]>, and include the main beam 


only. The values are taken from Planck Collaboration IV ( 2015) >, 
and are shown here to give an idea of typical beam shape at each 
frequency. Ellipticity is defined as the ratio of beam widths along 
the two main axes of the beam. A symmetric beam will have el¬ 
lipticity equal to 1. 
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Table 2. Table of beam characteristics. FWHM and ellip- 
ticity are based on the FEBeCoP effective beams ([Planck 
Collaboration IV 2015 1 . 3 earn efficiency includes the main and 


intermediate beam components. 


Frequency 

FWHM 

ellipticity 

efficiency 

30 GHz 

32.193” 

1.318 

0.99031 

44 GHz 

27.000” 

1.035 

0.99771 

70 GHz 

13.213” 

1.223 

0.99124 



Fig. 1. Squared symmetrized TT beam window for 30, 44, 70 
GHz. Solid: definition of Eq. ([7J. Dashed: definition of Eq. ([8]) 


For simulations we used the combined main and intermedi¬ 
ate radiometer scanning beams. The characteristics of these can 
be found in |Planck Collaboration IV| ( |2015| >. The last column of 
Table [2] gives the beam efficiency at each frequency. These are 
computed from of the harmonic beam expansion of the noise- 
weighted frequency beam. Beam efficiency is proportional to the 
monopole component of the expansion, and normalised so that 
the full An beam has efficiency 1. The squared value gives the 
effect in the spectral domain. The missing fraction reflects the 
power absorbed by far sidelobes. 

For comparison purposes we constructed a symmetrized 
beam window as follows. First we computed the weighted sum 
individual radiometer beams, 


bfm — 


z 



(6) 


where j labels the radiometers, and crj are their respective white 
noise rms. The m = 0 component of the beam expansion, 


Wf[m=0] 


47T 

2 { + \ 


?fl) 


(7) 


represents a beam that is symmetrized by averaging over the az¬ 
imuthal angle around the beam centre. Factor An/{21 +1) serves 
to normalize the window function to unity at ( - 0 for a beam 
with efficiency equal to 1. The power spectrum of the beam 
expansion provides an alternative way of constructing a sym¬ 
metrized beam, as proposed by Wu et alT[( )200l| >, 

Wqms] = 2C+\ 2 ^Cmf- (8) 



Fig. 2. TT spectrum at 30 GHz. We show the deconvolved spec¬ 
trum and map spectrum for white noise and realistic noise simu¬ 
lations. The spectra are averaged over 40 noise realizations to re¬ 
duce scatter. Unsmoothed deconvolved spectra rise towards high 
multipoles, roughly as proportional to the square of the inverse 
window function. The dashed line shows the uniform white noise 
spectrum divided by the squared beam window. 


where “ms” stands for mean of squares. The beam window of 
Eq. ([8]) is always larger than that from Eq. Q, thus representing 
a narrower beam. For LFI beams the difference is small. The 
squared window functions Wj for both definitions are plotted 
in Fig. [T] The squared window can directly be compared with 
power spectra. 

2.6. Results from frequency simulations 

Figure [2] shows the 30 GHz TT spectrum for white noise and 
realistic noise inputs. All spectra have been averaged over 40 
noise realizations. We show in the same plot the deconvolved 
spectrum and the binned map spectrum. The main effect from 
deconvolution is evident. While the white noise map spectrum is 
uniform, the deconvolved spectrum rises steeply towards higher 
multipoles. As a first approximation, the effect of deconvolution 
is to scale the binned map spectrum by the inverse of the sym¬ 
metric beam window squared (shown in dashed line type). This 
is, however, only true as a first approximation. The real deconvo¬ 
lution spectrum rises even more steeply due to beam asymmetry. 
We investigate this further in Sect. [3] 

The primary deconvolution output is the harmonic coeffi¬ 
cients axim- For many practical purposes one may want to trans¬ 
form the result into an ordinary sky map. The procedure of 
constructing a sky map from the deconvolved a X i m coefficients 
involves a smoothing operation that brings the noise spectrum 
back down to the level of un-deconvolved noise. We discuss this 
procedure in Sect. [6] In this and the following three sections, 
where we do analysis in harmonic space, we always show the 
unsmoothed raw noise. In harmonic space the smoothing is a 
simple (-dependent scaling operation which does not change the 
relative differences between spectra. The relative differences be¬ 
tween spectra are thus more important than their absolute level. 

Residual correlated noise dominates the lower multipoles, 
but decreases towards higher multipoles. Fig. [3] shows the ratio 
of correlated residual noise spectrum and white noise spectrum. 
The quantity shown is computed as (Cfr ea i - Cf wn )/Cf wn , where 
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Fig. 3. Ratio of the residual correlated noise spectrum and the 
white noise spectrum, for deconvolution (red) and ordinary map¬ 
making (blue). From top down: 30 GHz, 44 GHz, 70 GHz. 


Cfteai and C/ wn are the spectra from realistic and white noise sim¬ 
ulations, respectively. 

An interesting and perhaps surprising observation is that de- 
convolution leaves relatively less correlated residual noise on top 
of the white noise component at high multipoles, than does nor¬ 
mal mapmaking. This has important consequences. The level of 
correlated residual noise is a measure of the error we make if we 
neglect the correlated component when estimating the residual 
noise. In Sect. [3] we derive a model for the deconvolved white 
noise spectrum, which we use as an approximation for the full 
noise at high multipoles. The accuracy of the approximation can 
be read from Fig. [3] Further, MC simulations involving only 
white noise spectrum are much cheaper than simulations with 
realistic noise, as the destriping step is avoided. 

We now look closer at the white noise spectrum. The spectra 
vary strongly due to the limited number of realizations. However, 
when comparing spectra derived from the same simulation with 
different methods, we can make use of the fact that we have the 
same noise realizations in both cases. We take the white noise 
map spectrum as reference level, and divide the other spectra 
by it. This reduces the scatter in the spectra and shows the rela¬ 
tive differences clearly. Spectral ratios constructed this way are 
plotted in Figure [4] The white noise map spectrum is uniform 
by construction. The shape of the deconvolved white noise spec¬ 
trum reflects the beam window function. The deconvolved spec¬ 
trum rises above the binned map spectrum at the lowest mul¬ 
tipoles by approximately 2%. This is the effect of the missing 
power lost into the far sidelobe component of the beam. (The 
y =99% efficiency of Table [^translates to an 1 - y 2 w2% effect 
in spectral domain.) Deconvolution scales the signal up to cor¬ 
rect for the missing power, and the same increase is observed in 
the noise spectrum. 

The purple line on Fig. [4] is the result of a simulation where 
we deconvolved for a perfect delta beam (“deconvolution with¬ 
out deconvolution”). Even if the beam is a delta peak, the decon¬ 
volution procedure is not identical to the combined map-binning 
and anafast procedure. Deconvolution performs a maximum- 
likelihood fit of the harmonic coefficients to the TOI, while 
anafast computes a harmonic expansion over the sky map, 
without taking into account the varying hit counts. The two pro¬ 
cedures are identical only in the case where the distribution of 



1.05 


wn map 

— wn dec 

— delta 


0.95 


multipole 


30GHz TT 


Fig. 4. Ratio of the deconvolved 30 GHz white noise TT spec¬ 
trum, and the corresponding map spectrum. Orange: deconvo¬ 
lution with realistic beam. Purple: deconvolution with a delta 
beam. The deconvolved spectrum rises above the reference level 
due to missing beam power, which deconvolution compensates 
for. 

observations is uniform over the sphere. The spectrum from de- 
convolution with a delta beam is slightly below the reference 
level of the white noise map spectrum, but the difference is small. 
Above l = 600 the spectrum drops, which is a typical boundary 
effect in deconvolution. The drop occurs at a much higher mul¬ 
tipole than the drop in the spectral ratio of Fig. [3] The exact 
mechanism behind the latter drop remains unclear. 


2.7. Horn pair LFI18/23 

Finally we compare the two destriping options discussed in Sect. 


2.2 We performed two simulations, where we deconvolved the 


data from horn pair LFI18/23. In the first set of simulations we 
used for destriping the same 18/23 horn-pair data set, in the sec¬ 
ond one the full 70 GHz data set. The results are shown in in Fig. 
[5 We show the deconvolved spectra for white noise and for re¬ 
alistic noise, with both destriping options. The white noise map 
spectrum is shown in the same plot for reference. The scatter is 
large due to the small number of realizations (10), but the main 
conclusion is still clear. Full destriping leaves considerably less 
residual noise at low multipoles. 

In the lower panel we replot the spectra the same way we 
did for full 30 GHz data in Fig. [4] We divided the spectra by the 
white noise map spectrum to reduce scatter, and to bring out the 
differences more clearly. 


3. White noise covariance 

3. 1. White noise covariance 

We proceed to study the properties of deconvolved noise in more 
detail. The deconvolution operation of Eq. when operating 
on a noise time stream n, yields a vector of harmonic coefficients 
whose properties are described by the harmonic noise covari¬ 
ance matrix given by 

C = (A t C- 1 A)- 1 A t C; I <«ii r )C- I A(A t C; 1 A)- 1 . (9) 

It describes the noise correlation between elements a s i m . 


5 





















E. Keihanen et al.: Impact of beam deconvolution on noise properties 





(U 

(-0 

'o 


c 




Fig. 5. TT spectrum for horn pair LFI18/23. We show the spectra 
for devonvolved white and realistic noise. In the case of realistic 
noise we compare two destriping options: independent and full 
destriping. The plot is based on 10 Monte Carlo realizations. The 
lower panel shows the same spectra replotted in dimensionless 
units to bring out their relative differences; The undeconvolved 
white noise map spectrum (green) is taken as reference level, and 
all other spectra are divided by it. 


In this section we study the simpler case where the input TOI 
consists of pure white noise, with diagonal variance 



Fig. 7. Zoom to TT bias at 30 GHz and 70 GHz. Red: MC with 
realistic noise. Black: MC with white noise. Green (30 GHz): 
piecewise constructed noise bias. Dashed, estimate based on the 
beam window. The circles indicate the bias for ( =500 or 800, 
constructed from a partial NCV matrix of width (from down up) 
A£-5, 10, 20, 40 (30 GHz), or A£=5, 10, 20 (70 GHz) 



C n = (nn T ). (10) 

The harmonic noise covariance matrix simplifies into 

C = (A t C n 'Ar l . (11) 

The noise covariance is thus given by the inverse of the decon¬ 
volution matrix (A C,) 1 A). This is a product we can extract from 
artDeco. The size of the matrix rapidly increases with increas¬ 
ing ('max- The rank of the matrix with polarization is 3(Z max + l) 2 . 
The memory requirement increases as square of the rank, and the 
CPU time needed for inversion in the third power, or /^ lax . The 
full matrix inversion rapidly becomes prohibitive. In practice we 
can only construct it for a limited multipole range. 


Fig. 8. Effect of ( max parameter on 30 GHz TT white noise spec¬ 
trum. We deconvolve the same data set with f max =50, 100, 200, 
400. To reduce scatter, we divide each spectrum by the white 
noise map spectrum, and show them in dimensionless units. The 
plot is based on 40 Monte Carlo realizations. The smooth lines 
in same colour show the predicted noise bias, constructed from 
harmonic white noise covariance matrix, for (max =50, 100, 200. 
To bring them into units comparable with the Mont Carlo spec¬ 
tra, we divide them by the white noise level 3.2180 • 10~ 15 K 2 
(from Table [3j. The last 2-5 unreliable multipoles are excluded 
for clarity. 


3.2. Noise bias 

We want to compare the noise spectra from simulations to the 
covariance matrix. We derive from the covariance matrix a noise 
bias, a prediction for the expectation value of the noise spectrum 
of Eq. <(5j. 

We define gx s as factors that convert the spin harmonics a s ( m 
into a X f m . X =T,E, B. We write Eq. 0 as 


ClXCm 


^ ! gXs^stm 


( 12 ) 
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multipole i 


Fig. 6. TT, EE, and TE deconvolved noise spectrum for 30, 44, and 70 GHz. Black: white noise MC simulation. Red: realistic MC 
simulation. Dashed: estimate based on inverse symmetrized beam window. Green: noise bias from white noise covariance matrix. 
At 30 GHz we construct the bias for the complete multipole range, by a procedure that assembles the spectrum in pieces. At 44 and 
70 GHz we construct the bias for range C— 0-200, and individual values for multipoles £ = 300,400... 1400 


Now 


gT,0 = 1 

gE,2 = gE-2 = - 1/2 
gB2 = ~gB- 2 = 1 /2 
gT, 2 = gT-2 = gE, 0 = gB.O = 0. 


(13) 


The noise bias estimate is obtained from the noise covariance C 
as 

= Z (14) 


ss'mm' W 

Formally we can write 


where 


G 


XY( 

t' ms,l" m! s' 


C xn = Tr[CG™] 
1 


11 + 1 


gXsg Ys’$(('&({" Sri. 


(15) 

(16) 


The noise bias is extracted from a block-diagonal part of the 
NCVM, where each 3x3 block corresponds to a pair of (£, m) 
indeces, and ,v index takes all allowed values (0, +2). 

Figure [6] shows the TT, EE, and TE, noise spectra for three 
LFI frequencies. The BB spectra are quite similar to EE, and 
are not shown. We show the deconvolved white noise spectrum 
from MC simulation, along with a noise bias estimate obtained 
from the white noise covariance. We show also the realistic noise 
spectra, which rise above the white noise spectra at low multi¬ 
poles. The white noise component dominates the noise at high 
multipoles. Consequently, the white noise covariance provides a 
good model for the full noise at these high multipoles. We show 
also the raw estimate based on the symmetric beam window. It 
provides a reasonable approximation at intermediate multipoles, 
but underestimates the noise at high multipoles. 

Because of the large size of the deconvolution matrix, we 
cannot compute the exact noise bias for the full multipole range. 
We constructed the 30 GHz bias in pieces. First we constructed 
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and inverted the matrix for range £ = 0-110 to obtain the noise 
bias for multipoles £ - 0 - 100. We dropped the last 10 mul¬ 
tipoles to eliminate boundary effects. Similarly, noise bias for 
multipoles £ = 100 - 200 was constructed from a matrix com¬ 
puted for £ = 90 - 210. Multipoles £ = 200 - 300 were con¬ 
structed from two pieces in range £ = 190 - 260 and £ - 240 
- 310. Because the number of m multipoles increases with in¬ 
creasing £, we had to make the pieces narrower as we proceeded 
towards higher multipoles. Range £ = 300 - 400 was composed 
from pieces of 20 multipoles each, range £ - 400 - 600 from 
pieces of 10, and finally £ = 600 - 700 from pieces of 5 multi¬ 
poles. In all cases above £ - 300 we included 5 extra multiples at 
both ends, which we then dropped when combining the pieces. 
The whole computation took 20000 CPU hours. The estimate 
constructed this way is far more accurate than the one based on 
symmetric beam window, as can be seen from Fig. [6] However, 
the bias still underestimates the MC noise spectrum by 2 - 3% 
because it neglects correlations between distant multipoles. At 
30 GHz this is not a major problem, because we make an error 
of similar magnitude anyway by neglecting the correlated noise 
component (see Fig [3]). 

The computation becomes increasingly heavy towards 
higher multipoles. We therefore adopted a different procedure 
for 44 GHz and 70 GHz. We constructed the bias in one piece 
for multipoles £ - 0-200. At higher multipoles we concentrated 
the computational resources on a number of individual multi¬ 
poles, for which we estimated the bias as accurately as possible. 
We optimized the matrix inversion routine in such a way that 
it only computes the diagonal elements actually needed for the 
bias. We constructed the white noise covariance for the multi¬ 
poles from £() - Af to £o+A£, and picked the central value Iq for 
plotting. Values of £q were taken to be multiples of 100. 

At 44 GHz we used A£=30 for £ 0 = 300-500, and A£=20 for 
/o = 400 - 800. A fortunate coincidence comes to help here. The 
correlations between distant multipoles become weaker with in¬ 
creasing £, so that even though we had to make the pieces nar¬ 
rower, the error in the bias rather stayed constant. The noise bias 
values obtained for £q are shown by circles in Fig. [6] The dif¬ 
ferences between predicted noise bias and MC spectra are of the 
order of 1%, and in many cases the predicted value lies within 
the statistical scatter of the MC spectrum. 

At 70 GHz the computation is more demanding. We applied 
a similar procedure as at 44 GHz. We constructed the full co- 
variance for £ = 0 - 200, then evaluated the bias for multiples 
of 100 above that. We used A£=30 for Iq - 300-500, A£=20 for 
/o = 600 - 800, and A£=10 for Iq - 900 - 1400. We reach an 
accuracy of 1 - 3%. 

The importance of parameter A£ is illustrated in Fig. [7] We 
computed the noise bias for £ = 500 at 30 GHz, and for £ = 800 
at 70 GHz, for increasing values of A£. At 30 GHz we show 
the result for A£= 5, 10, 20, 40. The values approach the Monte 
Carlo result as A£ increases, as expected. The computational bur¬ 
den rapidly increases with increasing A£. While computing the 
A£=20 point took 94 minutes on 960 cores, the A£=40 point al¬ 
ready took three hours on 1920 cores. 

At 70 GHz we computed the bias for A£= 5, 10, 20. The last 
point took 4 hours on 1920 cores. Again the results are approach¬ 
ing the MC result, but the three data points are not sufficient to 
conclude if there is a real convergence. Also the statistical scat¬ 
ter in the MC spectrum, based on 10 realizations, is too large for 
us to give an accurate estimate of the remaining error, other than 
saying that it is of the order of 1%. 

The 1-3% error in bias estimation may still be too large for 
many applications. To compute the bias more accurately, we 


would need a more efficient matrix inversion routine, or one 
that constructs the diagonal of the inverse. Although more ef¬ 
ficient methods surely exist, they are not likely to bring us very 
far, since the largest matrices we invert already have rank over 
200 000. Should higher accuracy be required, the best approach 
is probably to produce a larger number of Monte Carlo simula¬ 
tions and to extract the bias from them. This will require signifi¬ 
cant optimization of the pipeline. 

The realistic noise simulations at 70 GHz took 15 min of 
wall-clock time per MC cycle on 960 cores, corresponding to 
250 CPU hours. Roughly half of that went into noise generation 
and destriping in Madam, another half into the deconvolution 
process. The cost of white noise simulations was 52 CPUh per 
realization. At 30 GHz, the cost was 52 CPUh for realistic and 
41 CPUh for white noise simulations. 

The limiting factor in our current pipeline is the disk space, 
rather than CPU time, since we are writing the intermediate 4D 
map products on disk. The pipeline could be made a lot more 
efficient by if the step of writing disk I/O could be avoided. In 
case of white noise simulations, this could be done by generating 
a noise 3D map directly in artDeco memory, the same way we 
generate a noise TOI in Madam. This remains a topic for future 
development. 

At 44 GHz, white noise simulations give sub-percent accu¬ 
racy above multipoles £ = 800, as can be seen from Fig. [3] At 
70 GHz, this accuracy is reached already around £ = 400. At 30 
GHz, the contribution from the correlated component is above 
2% over the whole multipole range. 

3.3. Resolution 

The auto spectra (TT, EE, BB) of a binned white noise map 
become uniform, when averaged over a large number of realiza¬ 
tions. The level of the spectrum is dependent on the hit count 
distribution. For a given total number of hits and fixed noise rms 
per sample, the lowest spectral level is obtained when the hits are 
uniformly distributed over the sky. Downgrading the resolution 
of a white noise map in general reduces the spectral level, as the 
distribution becomes more uniform. 

The spectral levels can be predicted from a known hit distri¬ 
bution, as described in |Planck Collaboration VI| ( j20151 >. We tab¬ 
ulate the predicted white noise levels for our simulation param¬ 
eters for different resolutions in Tabled We show the results for 
two radiometer weighting schemes: noise weighting and horn- 
uniform weighting. Horn-uniform weighting was used in actual 
LFI mapmaking. 

In deconvolution the situation is different. The resolution of 
the input 4D map does not have a significant effect on residual 
noise. Instead, resolution dependence enters through parameter 
(max- This is demonstrated in Fig. [8] We deconvolved the same 
30 GHz data set with £ max = 50, 100, 200, 400. We show also 
the corresponding noise bias estimates for £ max = 50,100, 200, 
constructed from the noise covariance matrix. For £ max =400 we 
were unable to compute the exact noise bias, due to the huge 
size of the covariance matrix. For plotting purposes the spec¬ 
tra have been divided by the spectrum of an A s id e =1024 (3.44’) 
white noise binned map. The noise bias are scaled to the same 
dimensionless units by dividing them by the theoretical white 
noise level of 3.2180 ■ 10 15 K 2 , taken from Table[3] The noise 
level at low multipoles increases with increasing £ max , as there 
are more multipoles to be solved from same amount of data. 

The result has relevance for the following section, where 
we construct the full low-resolution noise covariance matrix for 
multipole range £ = 0 - 50. In order to have an exact match be- 
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Table 3. White noise level for ordinary mapmaking, for noise 
weighting (nw) and horn-uniform weighting (huw), and for var¬ 
ious map resolutions (Wide)- The average size of a Wide =1024 
(64) pixel is 3.44’ (55’) 


C TT (10- 15 K 2 ) Cee.Cbb (10- 15 K 2 ) 


Aside 

nw 

huw 

nw 

huw 

30 GHz 

64 

3.1740 

3.1740 

6.3356 

6.3633 

128 

3.1910 

3.1914 

6.3666 

6.3944 

256 

3.2003 

3.2015 

6.3837 

6.4117 

512 

3.2067 

3.2083 

6.3971 

6.4251 

1024 

3.2180 

3.2199 

6.4294 

6.4576 

44 GHz 

64 

4.2407 

4.2917 

9.2814 

9.3898 

128 

4.2544 

4.3058 

9.3170 

9.4259 

256 

4.2620 

4.3135 

9.3364 

9.4455 

512 

4.2666 

4.3182 

9.3491 

9.4583 

1024 

4.2738 

4.3254 

9.3758 

9.4852 

70 GHz 

64 

3.3202 

3.3327 

6.6476 

6.7021 

128 

3.3291 

3.3432 

6.6626 

6.7172 

256 

3.3334 

3.3482 

6.6697 

6.7243 

512 

3.3360 

3.3511 

6.6745 

6.7292 

1024 

3.3381 

3.3535 

6.6808 

6.7355 


tween the covariance matrix and the actual a^fmcoefficients, it is 
required that both have been computed with same 7 rnax . 


4. Full noise covariance 


We now consider the more complicated case where the input 
consists of 1 // type noise, which has been destriped to remove 
the correlated component. We aim at deriving a full noise co- 
variance matrix, which takes into account the correlated noise 
component. We assume a destriping process that involves a noise 
prior. One such implementation is the Madam code presented in 
Keihanen et al. (20101. The map binned from the destriped TOI 
approaches that of generalised least squared (GLS) mapmaking 
methods (Poutanen et aL] |20061 |Ashdown et al.||20O7a|b| |2009| > 
when the baseline length approaches one sample. 


4.1. Destriping 


Consider then the residual noise n in a TOI stream that has 
been processed according to Eq. p0| ). We aim at deriving a for¬ 
mula for (nn T ), and finally for the covariance of Eq. (|9]», where 
h is now the destriped noise stream. We assume that the original 
noise stream n' obeys the model ( fI7| , so that 

</iV T > = C n + FC b F T . (21) 

The baseline length, hidden inside matrix F, is a key param¬ 
eter in the calculation. The baseline length plays a double role in 
the noise covariance. On the one hand, it limits the frequencies 
that are taken into account by the covariance matrix. The fre¬ 
quencies above the inverse of the baseline length are neglected. 
If we choose too long a baseline length, we will underestimate 
the residual noise. For optimal results, we should thus choose a 
baseline length as short as possible. 

On the other hand, the baseline length in F represents the 
baseline length that was used in destriping. If we select a shorter 
baseline for covariance computations, we will not model the 
destriping process correctly. However, simulations with Planck 
LFI data (Planck Collaboration VI 2015| > show that the destrip¬ 
ing result converges already around a baseline length of one sec¬ 
ond. There is correlated noise at higher frequencies, but destrip¬ 
ing is unable to remove it, and making the baseline shorter does 
not change the results. We can turn this around and argue that we 
do not make a significant error if we assume a shorter baseline 
length than was actually used for destriping, as long as the actual 
baseline length was short enough that the optimal solution was 
already reached. 

All this put together, we argue that we get the optimal noise 
covariance when we make as short a baseline as possible. We 
bring this to the extreme and set the baseline length to one TOI 
sample. With this assumption, matrix F reduces into a unity ma¬ 
trix, F = I. The destriped TOI ( f20| > simplifies into 

y = y - (C n -'z + ChVc-'Zy = (c b c n ‘z + iy l y . (22) 


4.2. Covariance matrix 


We can rewrite Eq. (22 1 for the destriped noise component 

n = (C b C^Z + 1 rv. 


(23) 


We reformulate this with the aim of constructing the time do¬ 
main covariance (nn T ). Inserting Z we obtain 

/! = C n [C n + C b - C b C; 1 P(P T C- 1 Pr 1 P T r 1 /i'. (24) 

We apply the Woodbury formula to get 


In the destriping approach the initial noise stream is modelled as 
ri = F b + w wn , (17) 


n = c n [c t - 1 + cr 1 c b c- 1 P[i-(P T c; 1 pr I p T c t - 1 c b c- 1 P]- 1 

x(P T C- l P)- l P T C; l ]n' (25) 


where « wn represents white noise, and the correlated component 
is modelled by a sequence of baselines b. Matrix F formally 
spreads the noise baselines into TOI. 

The destriping solution gives an estimate for the baseline 
vector b as 

b = (F r C-*ZF + C b l r l f T C- l Zy, (18) 

where 

z = i-c i ; 1 P(P r c [ ; 1 Pr 1 P r c n 1 . (19) 

Here P is the pointing matrix, and C b is the a priori covariance of 
the baselines, constructed from the known noise spectrum. The 
cleaned TOI is constructed as 

y= y -Fb=y- F(F r C-‘ZF + C^r'F^'Zy. (20) 


where 


C t - C n + C b 


(26) 


denotes the full (white+correlated) time domain noise covari¬ 
ance of un-destriped TOI. We bring (P T C~ 1 P) 1 inside the 
brackets and note that Cf'CbC” 1 = C n 1 - C, 1 . Equation (25 1 
simplifies to 


n = CnC.-V + C b C t - 1 P[P T C t - 1 Pr 1 P T C-V . 


(27) 


Assuming that the true noise obeys the noise model, we have 
C t = (n'n' T ), and we obtain for the time domain covariance of 
the destriped noise 

(nn T ) = G n G^G n + P[P T C^‘P]~‘P T (28) 

-C n C7 1 P[P T C t - 1 Pr 1 P T C-'C n . 
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Inserting this into Eq. (|9| we finally obtain for the harmonic 
noise covariance matrix 


C = (A t C“ 1 A) _1 A t C^ 1 A(A t C„ I A) _1 

+(A t C n 1 Ar 1 A t C n 1 P(P T C t ~ 1 P)~ 1 P T C l 7 1 A(A t C n 1 A)~ 1 
-(A t C; 1 Ar 1 A t C“ 1 P(P T Cr 1 P)- I P T C t r 1 A(A t C; I Ar I . (29) 


Equation (291 is the main result of this section. 


4.3. Cross-correlation 

As discussed in Sect. |2.2| there are two destriping options avail¬ 
able when deconvolving data subsets. Either we destripe the full 
data set (full destriping) or only the deconvolved data set (inde¬ 
pendent destriping). In the former case, the common destriping 
step generates correlation in residual noise between data sets. 
This effect too we can assess through the noise covariance for¬ 
malism. We can construct a cross-covariance matrix 

Ci2 = («1«2> (30) 

where ci\ and «b are the harmonic coefficients obtained by de- 
convolution of the two data sets. The cross-covariance matrix is 
constructed in very much the same way as the covariance matrix 
of Eq. ( [29] ), only now we have on both sides two destriping op¬ 
erations Ai and A 2 . Formally both operate on the same full data 
stream, but yield zeros when operating to a part of TOI that does 
not belong to the data set. If the two data sets do not overlap, we 
have 

A|C t " 1 A 2 = 0. (31) 

With this assumption, the cross-covariance matrix is given by 

C 12 - (AlC- 1 A 1 r 1 [AlC- I P(P T C t - 1 P)- I P T C; 1 A 2 

-A|C 7 1 P(P T C t - 1 Pr 1 P T C- 1 A 2 ](AjC- 1 A 2 r I . (32) 

The auto-covariance for the full destriping case is given by the 
same formula ( [29] ) as for independent destriping, only P is taken 
to be the pointing matrix of the full data set, while in case of in¬ 
dependent destriping it is interpreted as the submatrix that con¬ 
tains the rows for the data set in question. 


4.4. Implementation 

Equations ( [29] ) and ( [32] appear complicated, but actually consist 
of only five ingredient matrices that appear in different combina¬ 
tions: AtC-'A, A'C-'A, A'C^P, A'C-'P, and P'C^P. We can 
further reduce the computation load as we describe in the fol¬ 
lowing. As described in Keihanen & kcinecke] (2012j ), artDeco 
performs deconvolution through a process where the TOI is 
first compressed into a “3D map”. The data is binned on a 3- 
dimensional grid, spanned by the pointing angles In 6 

and <p, artDeco are pixelized as in HEALPix, and 1 jj is divided 
uniformly into A ps j cells. We can formally write the deconvolu¬ 
tion operation AC n 1 as 

A f C;' = GS T C n ', (33) 

where S T , operating on a TOI, represents the operation of bin¬ 
ning the TOI into a 3D map, and G the remaining part of the 
deconvolution operation. 

Matrix P T represents the operation of binning a TOI into a 
two-dimensional sky map. This too can be split into a two-step 
operation where one first constructs a 3D map, and from that a 
sky map. We write 

P T = X T S T . (34) 


An element of a 3D map is identified through a combined index 
qi, where q refers to sky pixel, and i = 0... /V psl - 1 to the bin 
defined by beam rotation angle t/r,. The elements of a sky map 
are referred to as Ip, Qp, or Up. In this notation, the X matrix is 
given by 

X qi.lp — &q£p 

~ COS (2((//, + lApol)) fiqep (35) 

X qi,Up — Si H ( 2((//, + (Apol )) tiqtp 

where <// po | is the detector-specific angle of polarization sensi¬ 
tivity. To add to the generality, we allow the 3D map pixel res¬ 
olution to be higher than that of the sky map. In this notation 
d qcp = if pixel q belongs into the larger pixel p. However, in our 
computations we always use same resolution for both, so that 6 
becomes the usual Kronecker delta. Operation X T is thus a triv¬ 
ial coaddition operation involving weighting by the polarization 
angle. 

The five ingredient matrices become 

A'C^A = G t (S T C“ 1 S)G (36) 

A'C^'A = G t (S T C7'S)G 

A t C‘ 1 P = G t (S T C‘ l S)X 

A t C‘ 1 P = G t (S T C t -'S)X 

P T C t 'P = X T (S T Cr‘S)X. 

In the end, there remain three matrices to evaluate: (S T C n 'S), 

(S T C t ( I S), and G. Operation X is carried out algorithmically 
where needed. Matrix G is part of the artDeco deconvolver, 
and we can extract it from the code. Matrix (S T C~' S) represents 
the white noise covariance between elements of a 3D map. It is 
diagonal, with elements given by cr 2 /N P i, where <x 2 is the white 
noise variance in TOI domain, and N p i is the number of hits to 
bin pi. 

To proceed with matrix (S T C t r 1 S), we first use the Woodbury 
formula to rewrite C" 1 as 

cr 1 = (Cn + Cb)- 1 = C- 1 - C-kCb 1 + C-'r'C- 1 . (37) 


This formulation has the advantage that C n appears as inverse, 
which makes dealing with flagged data easier. The operation of 
Eq. (37 1 can be interpreted as a filtering operation in time do¬ 
main. To construct one column (S T C t r l S), we pick a TOI from a 
3D map with one non-zero element, apply the filter of Eq. © 
to it, and compress the TOI back into a 3D map. The FFT tech¬ 
nique is usually efficient in filtering, but since we are dealing 
with TOIs with very few non-zero elements, we find more effi¬ 
cient a procedure where we first use FFT to construct the filter, 
but then to perform the actual filtering in time domain. When 
dealing with Planck data we assume that noise is uncorrelated 
from one pointing period (typically one hour) to another, and do 
the filtering by pointing period. 

We further benefit from the fact that matrices of Eq. (36 1 are 
additive. If we have two sections of TOI, with no noise correla¬ 
tion between, we can compute the contribution from each section 
to the ingredient matrices independently, and coadd the results. 
In particular, several detectors are combined simply by coadding 
their matrices. 

We construct the ingredient matrices for each radiometer 
quadruplet and each half-year survey, write them on disk, and 
then coadd them in different combinations in a separate step. 
That way we avoid doing repeated computations when dealing 
with different data combinations. 
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4.5. Flagged data 

To add to the generality of the method, we allow for a case where 
part of the data is flagged as unusable. Madam handles flags by 
setting formally C n 1 = 0 for the flagged samples. This is easy 
to take into account in (S T C n 'S). Matrix (S T C t r l S) requires one 
more manipulation step. Eq. ( |37} is expensive to construct, since 
the middle matrix C b 1 + C n 1 is no more stationary. We therefore 
make the approximation 

(C, 1 +C- 1 )- 1 «(C b ' +c;,')-' (38) 

where C wn is the original white noise variance, where flags have 
not been applied. This can be evaluated efficiently using Fourier 
techniques. We thus have for the last ingredient matrix 

S T C t -‘S * S T C n 'S - S'C/'/C, 1 + C^r'C-'S. (39) 

In absence of flags the formula is exact. The approximation mis¬ 
estimates the correlation between distant 3D map elements, but 
captures the main effect of flagging: the reduction in the amount 
of available data and the consequent increase in residual noise. 

4.6. Destriping mask and detector weighting 

The mapmaking procedure of Planck LFI involves a destriping 
mask ( [Planck Collaboration Vl|2015j ). The mask is used for the 
purpose of reducing errors arising from high signal gradients and 
from bandpass mismatch between detectors. 

The proper way of treating the mask in the NCVM would 
be to introduce another white noise covariance C lnas k. whose in¬ 
verse C“! k = 0 for samples that fall in the masked region. The 
masked covariance is then applied in the destriping phase, while 
the normal white noise covariance is applied in the deconvolu¬ 
tion phase. Unfortunately, unlike in the case of flags, the differ¬ 
ent covariance terms appear in combinations that do not cancel 
out, with the result that the noise covariance matrix of Eq. (|29| 
is no more valid, and not easily enhanced. Note that replacing 
the pointing matrix P by its masked version does not have the 
desired effect, since reducing P has the effect of reducing the 
residual noise, as there is a smaller number of pixels to solve, 
while the effect of destriping mask is to increase the residual 
noise level. 

The mask thus affects the residual noise in a way that is not 
captured by our noise covariance matrix. This has to be accepted 
as a source of uncertainty in noise modelling. 

In principle, a similar limitation arises from detector weight¬ 
ing. The destriping procedure of Planck FFI replaces the white 
noise covariance C„ by a slightly different modified covariance, 
which is made equal between a pair of detectors. This is referred 
to as horn-uniform weighting. The effect from this is expected to 
be small. 

4.7. Inputs 

The required inputs for NCVM computation are detector point¬ 
ing, harmonic beam expansion, and noise spectrum. It is not 
required that noise parameters remain the same throughout the 
mission, although in our simulations we have assumed so. The 
TOI can be assembled of sections with different noise spectra, 
and each section’s contribution added to the ingredient matrices. 

Four user parameters control the computation. Parameters 
(max and k max control the deconvolution process and enter matrix 
G. Parameters A 34 ? and /V 3D set the resolution of the 3D map di- 

side psi r 

mension of matrices S and G. The meaning of /V® e is twofold. 


Table 4. Fist of data combinations for which we construct the 
NCVM products. We give the data range ( Survey ) and combina¬ 
tion of horns (Horns) for the destriping phase and for the devon- 
volution phase. All products have been computed with decon¬ 
volution parameters f’ ma x=50, k max =2, N® e =64, W®=256. For 
each case listed we produce a harmonic noise covariance matrix, 
a N s ide=16 (220’) pixel-pixel covariance matrix, and a noise bias 
estimate. 


Deconvolution Destriping 


Input type 

Horns 

Survey 

Horns 

Survey 

30 GHz 

white 

27-28 

full 

27-28 

full 

wn+l/f 

27-28 

full 

27-28 

full 

44 GHz 

white 

24-26 

full 

- 


wn+l/f 

24-26 

full 

24-26 

full 

70 GHz 

white 

18-23 

full 

- 

- 

wn+l/f 

18-23 

full 

18-23 

full 

Single years 

white 

18-23 

yrl 

- 

- 

white 

18-23 

yr2 

- 


white 

18-23 

yr3 

- 

- 

white 

18-23 

yr4 

- 

- 

wn+l/f 

18-23 

yrl 

18-23 

full 

wn+l/f 

18-23 

yr2 

18-23 

full 

wn+l/f 

18-23 

yr3 

18-23 

full 

wn+l/f 

18-23 

yr4 

18-23 

full 

wn+l/f 

18-23 

yrl 

18-23 

yrl 

wn+l/f 

18-23 

yr2 

18-23 

yr2 

wn+l/f 

18-23 

yr3 

18-23 

yr3 

wn+l/f 

18-23 

yr4 

18-23 

yr4 

Horn pairs 

white 

18/23 

full 

- 

- 

white 

19/22 

full 

- 

- 

white 

20/21 

full 

- 

- 

wn+l/f 

18/23 

full 

18-23 

full 

wn+l/f 

19/22 

full 

18-23 

full 

wn+l/f 

20/21 

full 

18-23 

full 

wn+l/f 

18/23 

full 

18/23 

full 

wn+l/f 

19/22 

full 

19/22 

full 

wn+l/f 

20/21 

full 

20/21 

full 


On the one hand it sets the internal resolution of the deconvolu¬ 
tion calculation. As a rule of thumb, 1.5 should exceed f lnax . 
Thus, for deconvolution with £ max =50 one should have at least 
iV 3 ^ e =32, and for f rnax = 100 at least /V® e =64. On the other hand, 
V 3 ° e also sets the destriping resolution. The covariance matrix 
models a destriping procedure where the sky signal is assumed 
constant within a pixel of resolution /V S id e =(V® e . 

The required computational resources increase rapidly with 
increasing f max and /V 3[ <k . The complexity of the computation in¬ 
creases as Anax 5 , and the size of the matrix as (’max 4 . In practice 
we can produce a full covariance matrix for only the low multi¬ 
poles £ ^ 100. 

4.8. Constructing NCVM for Planck LFI frequencies 

We have computed NCVM products for Planck FFI frequencies, 
for a number of radiometer and survey combinations. A com¬ 
plete list is given in Table [4] For each radiometer and survey 
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Fig. 10. Structure of the TT block of the 30 GHz NCVM at low 
multipoles (£ = 0 - 6). Index on x and y axis is l 2 + £+m. Positive 
correlation is shown in red colour, negative in blue. 


combination we produced two sets of NCVM products: one with 
the white noise matrix, given by Eq. (Ill, and one with the full 
NCVM matrix given by Eq. 


__[l. We used the same pointing, 

noise parameters, and beam, as in the simulations of Sect. [2] 

We constructed the harmonic NCVM matrix for multipoles 
up to ^ max =50. Beam elements with high k only affect the high 
£ multipoles ( |Keihanen & Reinecke |2012| ). We chose k max =2 
for this low-multipole computation. Value f max =50 requires a 
minimum pixel resolution of V® e =32 (110’). To realistically 
model the mapmaking process of Planck LFI, we should have 
V® e =1024 (3.44’) , which is out of reach. As a trade-off, we 
used N® e =64 (55’) and (V®=256. The effect of this approxima¬ 
tion on matrix quality is assessed in Sect. [5] 

For all cases, the set of NCVM products includes: a harmonic 
NCVM matrix with £ max =50 (930 MB), a noise bias estimate 
computed according to Eq. and a pixel-pixel covariance 
matrix at resolution V s id e =16 (720 MB). All products include 
polarization. The construction of the pixel-pixel matrix from the 
harmonic matrix is described in Sect. [6] 

In addition to the full frequency matrices for 30, 44, and 70 
GHz, we produced matrices for 70 GHz horn pairs 18/23, 19/22, 
20/21, and for individual years (years 1-4) involving all 70 
GHz radiometers. For each of these, we considered both full and 
independent destriping. 

The construction of the ingredient matrices as described in 
Sect. 4.4 took 2000 CPU hours for each year and horn pair, 
adding up to 24000 CPU hours for the whole 70 GHz data set. 
From those we constructed the NCVM for different horn pair and 
year combinations. The process of constructing the final NCVM 
from the ingredients was dominated by the inversion of the mid¬ 
dle matrix, and took about 1500 CPU hours per matrix. 

The predicted noise bias for FFI frequencies is shown in Fig. 
[9] We show three auto spectra (TT, EE, BB) and three cross 
spectra (TE, TB, EB ). The cross spectra are small, and are mul¬ 
tiplied by 10 in the plot to show the structure more clearly. 

Figure [T0| illustrates the matrix structure in t, m space. We 
show the real part TT block of the 30 GHz NCVM for lowest 
multipoles (£ < 7). The matrix elements are arranged according 
to the combined index i = £ 2 + £ + in. The strongest correlations 
appear between elements with same £, or same m and £ separated 
by an multiple of 2. 


5. Validation of the covariance matrix 

5.1. Low-resolution simulations 

We validated our NCVM products through MC simulations. A 
complete list of simulations is given in Table [ 5 ] The simulations 
followed the same two-step procedure as the high-resolution 
simulations reported in Sect. [2] We first generated and destriped 
a noise stream with the Madam destriping code, and then ran the 
artDeco deconvolution code, to yield an array axim of harmonic 
coefficients. The noise and beam parameters are the same as in 
Sect. [2] This time we ran the deconvolution step with parame¬ 
ters f max —50, k max =2, in line with the settings used in NCVM 
construction. With this low an £ max value the deconvolution step 
is fast, and the computation time is dominated by generation of 
noise rather than deconvolution. We produced 100 noise realiza¬ 
tions for each case studied. As before, we considered two types 
of input noise: white and realistic 1 // noise. 

We cannot expect the NCVM to perfectly model the residual 
noise, due to simplifications we have to do when constructing 
the matrix, either due to limitations of the formalism, or due to 
limited computing resources. We can list four factors: 

1. The covariance matrix assumes destriping resolution of 
V s jde=64, while actual FFI mapmaking uses (V s id e =1024. 

2. FFI mapmaking involves a destriping mask, neglected by the 
NCVM. 

3. FFI mapmaking applies horn-uniform weighting in the de¬ 
striping phase, NCVM computation intrinsically uses noise 
weighting. 

4. NCVM assumes a baseline length of one sample, while the 
actual baseline length is 1.0 or 0.25 sec. 

We performed two series of simulations, which we refer to 
as “ideal” and “realistic”. The destriping parameters are given 
in Table [5] The realistic simulations mimick the FFI mapmak¬ 
ing procedure as closely as possible: we use a destriping mask, 
apply horn-uniform weighting, and destripe at a high resolution 
Wide = 1024. 

In the ideal simulations, the destriping parameters follow the 
intrinsic assumptions of NCVM. Destriping is performed at low 
resolution N s a e =64, no destriping mask is applied, and horn- 
uniform weighting is replaced by noise weighting. To exactly 
match the assumptions of NCVM computation, we ought to use 
a baseline length of one sample. This becomes, however, compu¬ 
tationally very demanding. As a trade-off we selected a baseline 
length of a quarter of the one used in actual mapmaking. 

The aim of the ideal simulation set is to verify that the 
NCVM is constructed correctly. The realistic simulation set tests 
how accurately the matrix produced models the actual FFI data 
products, given the approximations involved. 


5.2. Chi2 tests 


We validate the noise covariance matrices through^- 2 tests. We 
generated 100 realizations of of harmonic vectors axt m , as de¬ 
scribed above. For each vector and corresponding matrix we 
compute the value 


X 


2 


a T C l a 

Ndof 


(40) 


Here C is the covariance matrix, and Ndo f is the number of de¬ 
grees of freedom. Multipole £ contributes 'll + 1 degrees of free¬ 
dom, except for the monopole and dipole, which are identically 
zero for E and B. This yields in total Vd 0 f=3(f’ max +1) 2 - 8 de¬ 
grees of freedom, which for f max =50 gives Vdof=7795. For a 
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Table 5. List of MC simulation cases. The parameters from left to right are: Simulation identifier; input data type; combination 
of horns and survey included in deconvolution phase (“full” refers to the full 4-year data set, “yrN” to single years N=l,2,3,4); 
deconvolution parameters £ max and k max ; and N™ resolution of the 3D map; horn combination and survey in destriping phase; 
baseline length (s); destriping resolution; destriping mask used (T/F); detector weighting (nw=noise weighting, huw=hom-uniform 
weighting); number of Monte Carlo realizations. 


Deconvolution Destriping 


ID 

data 

Horns 

Surv. 

Anax 

^max 

AT® 

side 

N™ 

psi 

Horns 

Surv. 

Base 

^side 

Mask 

Weight 

MC 

30 GHz 

white 

white 

27-28 

full 

50 

2 

64 

256 

27-28 

full 

- 

- 

- 

- 

100 

ideal 

wn+l/f 

27-28 

full 

50 

2 

64 

256 

27-28 

full 

0.063 

64 

F 

nw 

100 

real 

wn+l/f 

27-28 

full 

50 

2 

64 

256 

27-28 

full 

0.25 

1024 

T 

huw 

100 

hres wn 

white 

27-28 

full 

700 

6 

1024 

4096 

27-28 

full 

- 

- 

- 

- 

40 

hres real 

wn+l/f 

27-28 

full 

700 

6 

1024 

4096 

27-28 

full 

0.25 

1024 

T 

huw 

40 

44 GHz 

white 

white 

24-26 

full 

50 

2 

64 

256 

24-26 

full 

- 

- 

- 

- 

100 

ideal 

wn+l/f 

24-26 

full 

50 

2 

64 

256 

24-26 

full 

0.25 

64 

F 

nw 

100 

real 

wn+l/f 

24-26 

full 

50 

2 

64 

256 

24-26 

full 

1.0 

1024 

T 

huw 

100 

hres wn 

white 

24-26 

full 

1000 

6 

1024 

4096 

24-26 

full 

- 

- 

- 

- 

10 

hres real 

wn+l/f 

24-26 

full 

1000 

6 

1024 

4096 

24-26 

full 

1.0 

1024 

T 

huw 

10 

70 GHz 

white 

white 

18-23 

full 

50 

2 

64 

256 

18-23 

full 

- 

- 

- 

- 

100 

ideal 

wn+l/f 

18-23 

full 

50 

2 

64 

256 

18-23 

full 

0.25 

64 

F 

nw 

100 

real 

wn+l/f 

18-23 

full 

50 

2 

64 

256 

18-23 

full 

1.0 

1024 

T 

huw 

100 

hres wn 

white 

18-23 

full 

1500 

6 

1024 

4096 

18-23 

full 

- 

- 

- 

- 

10 

hres real 

wn+l/f 

18-23 

full 

1500 

6 

1024 

4096 

18-23 

full 

1.0 

1024 

T 

huw 

10 

Single years 
white 

white 

18-23 

yrN 

50 

2 

64 

256 

18-23 

yrN 





100 

ideal indp 

wn+l/f 

18-23 

yrN 

50 

2 

64 

256 

18-23 

yrN 

0.25 

64 

F 

nw 

100 

real indp 

wn+l/f 

18-23 

yrN 

50 

2 

64 

256 

18-23 

yrN 

1.0 

1024 

T 

huw 

100 

ideal full 

wn+l/f 

18-23 

yrN 

50 

2 

64 

256 

18-23 

full 

0.25 

64 

F 

nw 

100 

real full 

wn+l/f 

18-23 

yrN 

50 

2 

64 

256 

18-23 

full 

1.0 

1024 

T 

huw 

100 

Horn pairs 
white 

white 

18/23 

full 

50 

2 

64 

256 

18-23 

full 





100 

white 

white 

19/22 

full 

50 

2 

64 

256 

18-23 

full 

- 

- 

- 

- 

100 

white 

white 

20/21 

full 

50 

2 

64 

256 

18-23 

full 

- 

- 

- 

- 

100 

ideal indp 

wn+l/f 

18/23 

full 

50 

2 

64 

256 

18/23 

full 

0.25 

64 

F 

nw 

100 

ideal indp 

wn+l/f 

19/22 

full 

50 

2 

64 

256 

19/22 

full 

0.25 

64 

F 

nw 

100 

ideal indp 

wn+l/f 

20/21 

full 

50 

2 

64 

256 

20/21 

full 

0.25 

64 

F 

nw 

100 

real indp 

wn+l/f 

18/23 

full 

50 

2 

64 

256 

18/23 

full 

1.0 

1024 

T 

huw 

100 

real indp 

wn+l/f 

19/22 

full 

50 

2 

64 

256 

19/22 

full 

1.0 

1024 

T 

huw 

100 

real indp 

wn+l/f 

20/21 

full 

50 

2 

64 

256 

20/21 

full 

1.0 

1024 

T 

huw 

100 

ideal full 

wn+l/f 

18/23 

full 

50 

2 

64 

256 

18-23 

full 

0.25 

64 

F 

nw 

100 

ideal full 

wn+l/f 

19/22 

full 

50 

2 

64 

256 

18-23 

full 

0.25 

64 

F 

nw 

100 

ideal full 

wn+l/f 

20/21 

full 

50 

2 

64 

256 

18-23 

full 

0.25 

64 

F 

nw 

100 

real full 

wn+l/f 

18/23 

full 

50 

2 

64 

256 

18-23 

full 

1.0 

1024 

T 

huw 

100 

real full 

wn+l/f 

19/22 

full 

50 

2 

64 

256 

18-23 

full 

1.0 

1024 

T 

huw 

100 

real full 

wn+l/f 

20/21 

full 

50 

2 

64 

256 

18-23 

full 

1.0 

1024 

T 

huw 

100 

hres wn full 

white 

18/23 

full 

1500 

6 

1024 

4096 

18-23 

full 

- 

- 

- 

- 

10 

hres real full 

wn+l/f 

18/23 

full 

1500 

6 

1024 

4096 

18-23 

full 

1.0 

1024 

T 

huw 

10 

hres wn indp 

white 

18/23 

full 

1500 

6 

1024 

4096 

18/23 

full 

- 

- 

- 

- 

10 

hres real indp 

wn+l/f 

18/23 

full 

1500 

6 

1024 

4096 

18/23 

full 

1.0 

1024 

T 

huw 

10 

hres wn full 

white 

19/22 

full 

1500 

6 

1024 

4096 

18-23 

full 

- 

- 

- 

- 

10 

hres real full 

wn+l/f 

19/22 

full 

1500 

6 

1024 

4096 

18-23 

full 

1.0 

1024 

T 

huw 

10 

hres wn indp 

white 

19/22 

full 

1500 

6 

1024 

4096 

19/22 

full 

- 

- 

- 

- 

10 

hres real indp 

wn+l/f 

19/22 

full 

1500 

6 

1024 

4096 

19/22 

full 

1.0 

1024 

T 

huw 

10 

Special simulations 

delta beam wn 

27-28 

full 

700 

2 

1024 

4096 

27-28 

full 





40 

delta beam 

wn+l/f 

27-28 

full 

700 

2 

1024 

4096 

27-28 

full 

0.25 

1024 

T 

huw 

40 

Anax 

wn 

27-28 

full 

100-400 

2 

1024 

4096 

27-28 

full 

- 

- 

- 

- 

40 

Anax 

wn+l/f 

27-28 

full 

100^100 

2 

1024 

4096 

27-28 

full 

0.25 

1024 

T 

huw 

40 

delta beam 

white 

18-23 

full 

1500 

2 

1024 

4096 

18-23 

full 

- 

- 

- 

- 

10 

delta beam 

wn+l/f 

18-23 

full 

1500 

2 

1024 

4096 

18-23 

full 

1.0 

1024 

T 

huw 

10 

Anax 

white 

18-23 

full 

100-800 

2 

1024 

4096 

18-23 

full 

- 

- 

- 

- 

10 

Anax 

wn+l/f 

18-23 

full 

100-800 

2 

1024 

4096 

18-23 

full 

1.0 

1024 

T 

huw 

10 

Ws.de 1024 

wn+l/f 

18, 23 

full 

50 

2 

64 

256 

18-23 

full 

0.25 

1024 

F 

nw 

50 

mask 

wn+l/f 

18, 23 

full 

50 

2 

64 

256 

18-23 

full 

0.25 

64 

T 

nw 

50 

weight 

wn+l/f 

18, 23 

full 

50 

2 

64 

256 

18-23 

full 

0.25 

64 

F 

huw 

50 

base 

wn+l/f 

18, 23 

full 

50 

2 

64 

256 

18-23 

full 

1.0 

64 

F 

nw 

50 
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Fig. 9. Predicted noise bias for LFI frequencies: TT, EE, BB, TE, TB, EB. The cross spectra have be multiplied by 10 to bring out 
the structure more clearly. 
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Fig. 11. Reduced^ 2 distribution for 30 GHz, 44 GHz, 70 GHz, 
for white noise (left), ideal, and realistic simulation (right), for 
100 realizations. Also shown is the^ 2 distribution corresponding 
to an ideal result. The numbers within the plot give the mean 
value of the distribution. 


matrix that perfectly describes the noise properties, the results 
obey the Chi distribution. For large values of Ndof the distribu¬ 
tion becomes Gaussian, with expectation value 1 and standard 
deviation of y/2/Nd 0 f, which for AW=7795 is 0.016. 

We compared the white noise covariance matrix of Eq. o 
against white noise simulations, and the full noise matrix of Eq. 
( [29] ) against ideal and realistic simulations. Results from tests at 
the three Planck LFI frequencies, 30, 44, and 70 GHz are shown 


in Fig. 11 We show a histogram of the x~ values for 100 real¬ 


izations, and quote their mean. The white noise matrix models 
the white noise residual nearly perfectly, x 2 mean ranging from 
0.9992 to 1.0004, Comparison between the full matrix and ideal 
simulation gives slightly larger values, from 1.0002 at 70 GHz 
to 1.0037 at 44 GHz, while the expected one-sigma deviation for 


100 realizations is 0.0016. Deviations from one are larger in the 
case of realistic simulation, as expected, given the approxima¬ 
tions involved. At 30 GHz we find a x 2 mean of 1.0196, at 70 
GHz 1.0037. 

Results for 70 GHz subsets of are shown in Fig. [12] In the 
first column we compare the white noise matrix against a pure 
white noise simulation. Second and third column show results 
for the full matrix. Results are shown for two destriping options: 
full 70 GHz destriping, and independent destriping. 

In all cases, we find very good agreement with the white 
noise matrix and corresponding simulations. The same can be 
said about the full matrix and the ideal simulations. The^f 2 mean 
is at maximum 1.0038. In the case of realistic simulation, devi¬ 
ations are again larger. The agreement is, however, still quite 
good for cases where deconvolution and destriping involve the 
same data (independent destriping), the maximal x 2 value being 
1.0079. 

The case of full destriping, where the whole 70 GHz data 
set is destriped together, but only a subset of the cleaned data 
is provided to deconvolution, turns out to be more difficult to 
model with the covariance matrix. While the ideal simulation 
still shows good agreement, the mean \ ,2 results for the realis¬ 
tic simulation range from 1.0146 to 1.0495. The discrepancy is 
largest for horn pair LFI18/23, which we now pick up for a closer 
examination. 


We performed yet another series of simulations, where we 
applied full destriping, and deconvolved the LFI 18/23 horn pair. 
We took as starting point the ideal simulation, and turned on, 
one at a time, each of the features that distinguish the realis¬ 
tic simulation from the ideal one. That is: destriping resolution 
jVside=1024 instead of 64, horn-uniform weighting instead of 
noise weighting, baseline length 1 s instead of 0.25 s, and de¬ 
striping mask. We generated 50 noise realizations for each case. 
The results of the^ 2 test for these simulations are shown in Fig. 


13 We observe that detector weighting and destriping resolution 


have very little effect on the result. The finite baseline length and 
destriping mask both contribute significantly. Both effects alone 
raise the^ 2 mean to 1.03. 


As discussed in Sect. |4.6[ the mask affects the residual noise 
in a way that we cannot model in the covariance matrix. The ef¬ 
fect of baseline length, instead, could be cured by running map¬ 
making with a shorter baseline. 
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Fig. 12. Reduced^ 2 distribution for 70 GHz subsets. From top to 
bottom: full (4 years) 70 GHz data, years 1-4 (full frequency), 
horn pairs 18/23, 19/22, 20/21 (4 years). From left to right: white 
noise only, independent destriping ideal/realistic, full destriping 
ideal/realistic. In the case of full mission, independent and full 
destriping options are identical. 



Reduced x 2 


Fig. 13. Effect of various approximations on the^ 2 result for the 
worst LFI18/23 case. We add to the ideal simulation one non¬ 
ideality at a time: destriping mask, horn-uniform weighting, fi¬ 
nite baseline length, and destriping resolution. The result for the 
ideal simulation case is shown in the background. 



□ t =100 


I—I =200 




1.00 1.05 
Reduced x 2 
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Fig. 14. Effect of f m . IX parameter on x 1 test. We deconvolve the 
full 70 GHz data set (10 realizations) with (max =50, 100, 200, 
400, 800, 1500, and compare the lowest t = 0-50 multipoles 
with the NCVM with (max =50. 


5.3. Effect of ( max 

The noise covariance matrix models a deconvolution process 
where the harmonic coefficients are solved up to multipole ( lnax . 
For a good match between the matrix and the harmonic coef¬ 
ficients, the matrix and the actual deconvolution must use the 
same value of / max . One might prefer running deconvolution 
with a high £ mm to avoid aliasing effects, and then pick the lowest 
multipoles for a low-resolution data set. This, however, degrades 
the match between the matrix and the data. 

We illustrate this in figure [14] We deconvolved the 10 re¬ 
alizations of our 70 GHz full noise simulation with different 
(maxi pick the lowest ( = 0-50 multipoles, and ran a x 2 
test with the 7’ m a X =50 NCVM. The matrix gives a good match 
with the coefficients obtained with £ max -50, but at (’ m ax=100 
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Fig. 16. Noise bias for 30 and 44 GHz. TT and EE spectra com¬ 
pared against MC simulations. Line types are the same as in Fig 


E3- 



Multipole i 


Fig. 17. Noise bias in TT and EE for horn pair LFI18/23 com¬ 
pared against MC simulations: Solid: full destriping. Dashed: 
independent destriping 


there is already a significant deviation. Going from f’ max =200 
to f’ max =1500 changes the result only little. This is in line with 
results obtained for white noise, shown in Fig. [8] 

It is beyond the scope of this work to answer the question 
if the low fmax can cause signal distortion. Should this be the 
case, we still have a few options left. We can improve the match 
by constructing the NCVM to higher Computational cost 
increases rapidly with ^ max , but fm ax =50 used in this work is by 
no means an upper limit. Luckily, figure [14] suggests that the 
resolution effect has already largely saturated at ( max = 100, so 
increasing f max even a little may already solve the problem. 

Another option worth exploring is to construct the full matrix 
as combination of two components, a white noise component, 
computed to higher resolution (as we have done in Sect. [3]» and 
from a full noise component computed to lower resolution. 


Figure [TT] shows a comparison for horn pair LFI18/23. We 
consider again two destriping options. Simulations show that full 
destriping yields a much lower residual noise level, and this is 
captured by the NCVM as well. 

Figure |T8| shows the cross spectrum between 70 GHz horn 
pairs 18/23 and 19/22. When the horn pairs are destriped inde¬ 
pendently, there is no correlation in the noise. When the whole 
70 GHz data set is destriped together, there is a significant cross¬ 
correlation at low multipoles. The NCVM models this correctly. 
We show in the same plot also the auto spectrum for horn pair 
18/23 (same as in Fig. [TT}. The theoretical noise bias agrees 
well with the idealised simulation, but underestimates the real¬ 
istic noise by a few per cent, in agreement with the x 2 result for 
the same case. 


5.4. Validation through noise bias 

Another way of validating the NCVM product is to compare the 
predicted noise bias of Eq. ( fl5| against the MC spectrum of Eq. 
([5]). Figure [15] shows the comparison for the full 70 GHz data. 
We show all six spectral components. The MC spectra shown 
are averaged over 100 noise realizations. We show both the ide¬ 
alised and the realistic simulation. For auto-spectra (TT, EE, 
BB ) we show also the white noise simulation and the corre¬ 
sponding noise bias. 

The predicted auto-spectra agree reasonably well with the 
simulation, except for multipole 1 = 2, where the noise bias un¬ 
derestimates the true noise level. In the case of the much smaller 
cross spectra (TE, TB, EB), the 100 realizations do not provide 
enough statistics for a meaningful comparison. In the subsequent 
analysis we therefore concentrate on TT and EE spectra. 

In Fig. [16] we show the TT and EE spectra for 30 and 44 
GHz channels. The agreement is again good, and the deviation at 
( = 2 is smaller than in case of 70 GHz. The difference between 
ideal and realistic simulations is small. 


6. Covariance in pixel domain 

The primary output of artDeco is an array of harmonic coef¬ 
ficients. Consequently, the primary noise covariance matrix is 
constructed in harmonic space. 

From the harmonic coefficients one can further construct a 
sky map via spherical harmonic transform. The harmonic coeffi¬ 
cients must be multiplied by a smoothing kernel before perform¬ 
ing the expansion. Smoothing serves two purposes. Firstly, it re¬ 
moves ringing artefacts that arise around point sources and other 
regions of high signal gradients. Secondly, smoothing compen¬ 
sates for the rise in noise level towards high multipoles. 

6.1. High-resolution maps 

We consider first high-resolution maps. We do not attempt to 
build a covariance matrix for high resolution, but study the noise 
properties through Monte Carlo simulations. 

We transformed our high-resolution simulations to pixel 
space. Eliminating ringing artefacts in general requires a 
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Fig. 15. 70 GHz noise bias compared against MC simulations. We compare the noise bias from the full NCVM (orange) against ideal 
(light blue) and realistic (dark blue) simulations. For the auto spectra we show also the white noise bias (red) and corresponding 
MC simulations (green). The MC spectra are averaged over 100 noise realizations. 
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Fig. 18. TT spectrum for 18/23 and 19/23 (70 GHz) data sets. 
Each group of three lines shows a comparison between MC sim¬ 
ulations and corresponding NCVM prediction. The NCVM pre¬ 
diction is shown in orange or brown, the MC results in blueish 
colours. The topmost group shows the auto spectrum of 18/23 
(blow-up of Fig. [TT] ). The group in the middle is the cross¬ 
spectrum between 18/23 and 19/22, for full destriping. The 
group shown in dashed line type shows the cross-spectra for in¬ 
dependent destriping. 



Fig. 19. Deconvolved TT spectrum after smoothing with a 
Gaussian beam of FWHM=40’ (30 GHz), 30’ (44 GHz) or 20’ 
(70 GHz). Smoothing suppresses the noise at high multipoles. 
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Fig. 20. A 1 Jy point source and noise at same location 
(0.0373°,0) in a simulated 30 GHz map. We show a patch of 2° 
around the source. Top: binned map. Bottom: deconvolved and 
smoothed to 40’. Deconvolution stretches noise in direction per¬ 
pendicular to beam elongation. Smoothing suppresses the noise 
amplitude. 


smoothing width larger than the average width of the actual 
beam. The same conclusion is supported by noise analysis. As 
shown in Sect. [3j the deconvolved noise spectmm rises more 
steeply than the inverse symmetrized beam window. Smoothing 
by this window will thus not be sufficient. 

We applied a Gaussian smoothing kernel of width 
FWHM=40’ (30 GHz), 30’ (44 GHz), or 20’ (70 GHz). The 
widths were chosen so as to remove ringing effects at each re¬ 
spective frequency. Figure p~9] shows the TT noise spectra with 
smoothing applied. All spectra vanish at high multipoles. 

Figure [20] illustrates the effect of deconvolution on noise in 
the pixel domain. We show a simulated 1 Jy point source at 
30 GHz, located just above the origin of the Galactic coordi¬ 
nate system. In order to have the source in the centre of an 
M;ide=1024 pixel we place it not at the exact Galactic centre, 
but at (0.0373°,0) (latjon), which is the centre of the pixel just 
above. We show also as example one realization of noise at the 
same location. The noise in the binned map is dominated by 
white noise. The smoothing operation associated with deconvo¬ 
lution suppresses the noise level in the deconvolved map. 

Owing to the smoothing, the image of a point source on a 
deconvolved map is not point-like. It is, however, perfectly sym¬ 
metric. At the same time the deconvolution procedure generates 
correlations in the noise component, in the direction perpendic¬ 
ular to the direction of beam main axis. To quantify this, we 
compute an estimate of the correlation of the central pixel with 
its surroundings, as 


CmapXYj - (nxiV-Yj) (41) 

where X, Y = /, Q,U are the Stokes components,;; ; is the noise 
in pixel index i represents the reference pixel, and j labels all 



—60 0 60 

(^K ) 2 

Fig. 21. Correlation structure between the central pixel and its 
surroundings, computed from 40 realizations. We show the same 
sky region as in Fig. 20 



Fig. 22. 70 GHz II correlation between a reference pixel located 
at (2.388°,0) and the rest of the sky. 


other pixels of the sky. We evaluate C map as an average over the 
40 available noise realizations. The II component map is plotted 
in Fig. 22 As expected, noise is correlated in a direction at right 
angles to the beam main axis. 


6.2. Low-resolution maps and pixel NCVM 


We proceed to produce low-resolution maps, and related covari¬ 
ance matrices. We transformed the low-resolution simulations of 
Sect.[5]to pixel space. We adopted the smoothing procedure that 
was applied to LFI 2015 release data (Planck Collaboration VI 
2015). We applied Gaussian smoothing with FWHM=440’ and 
constructed maps at resolution A' SK | L ,= 16 (220’). 

The same operations must be applied to the covariance ma¬ 
trix. The smoothing operation reduces the effective number of 
degrees of freedom, making the covariance matrix singular or 
badly scaled. This is cured by adding a small amount of white 
regularization noise within each pixel of the map. A correspond¬ 
ing constant variance was added on the diagonal of the covari¬ 
ance matrix. Following the LFI procedure, we added white reg¬ 
ularization noise with rms=2 pK (I) or rms=0.02 pK (Q and U ) 
on each pixel. 

The operation on map can formally be written as 


m = Ha + « reg (42) 

Here H represents the harmonic expansion combined with 
smoothing, a is the vector of harmonic coefficients, as produced 
by artDeco, m is the beam-deconvolved map, and n reg repre¬ 
sents regularization noise. 
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The noise in map of Eq. ( |42] > is characterized by the covari¬ 
ance matrix 

C pix = HCH t + N reg (43) 

where C is the harmonic covariance matrix, and N leg is diagonal 
matrix representing the regularization noise. 

The operation HCH T is computationally intensive and bro¬ 
ken down into the following individual steps: 


1. Each row of C is transformed via SHT. This is done using the 
efficient lib sharp implementation (Reinecke & Seljebotn 
|2013) >. Since this library only supports real-valued maps, two 
SHTs (for the real and imaginary part, respectively) must be 
carried out for each row; depending on the parameter choice, 
these are either scalar or fully polarized. Since a matrix row 
is located entirely on a single MPI task, no communication 
is necessary for this step. 

2. To perform the analogous operation on every column of 
the intermediate result C H J , this matrix must be effectively 
transposed in memory, so that columns reside on a single 
MPI task., as the rows did before. This is a communication¬ 
intensive step. 

3. SHTs are carried out over the rows of (C H T ) T in a fashion 
analogous to step 1, resulting in C 7 x . 

4. In another communication step, this matrix is again trans¬ 
posed. 


Figure [22] illustrates the structure of the pixel covariance 
matrix. We show as a map one column of the 70 GHz covari¬ 
ance matrix, which represents II correlation between the refer¬ 
ence pixel and the rest of the sky. We pick the reference pixel 
just above the galactic centre. The reference pixel is centred 
at (2.388°,0) and has pixel number 1440 in the ring pixeling 
scheme of Healpix at resolution Wide=16- 

The reference pixel is strongly correlated with the neigh¬ 
bouring pixels. A weaker correlation stmcture follows the scan¬ 
ning direction. The map has a significant positive minimum of 
roughly 0.012 /./K 2 , reflecting the difficulty of distinguishing the 
map monopole from a global noise offset. The correlation struc¬ 
ture may vary a lot according to the reference pixel chosen. The 
plot shown here represents one example. 


6.3. x 2 tests in pixel space 

We validate the pixel covariance matrices through x" 2 tests, sim¬ 
ilarly to the validation tests presented in Sect. [5] The number of 
degrees of freedom is now given by the number of pixels in the 
Healpix map, times 3 Stokes components, which for Wide = 16 
gives /V dl) | =9216. Consequently, the one-sigma statistical varia¬ 
tion becomes 0.0147, and rms of mean for 100 realizations be¬ 
comes 0.00147. 

As in Sect. [5] we compared the white noise covariance ma¬ 
trix against white noise simulations, and full covariance matrix 
against ideal and realistic simulations. The x 2 results for full 
four-year frequency maps and for 70 GHz subsets are shown in 
Figs. [23] and [34] 

The mean reduced^ 2 values for pixel covariance are in gen¬ 
eral closer to one than in the case of harmonic covariance. This 
is the effect of regularization noise. We are adding to the maps a 
noise component that is accurately modelled by the extra term in 
covariance matrix, thus improving the agreement with the matrix 
and the maps. 

The x 2 mean is 1.01 for 30 GHz 4-year data, 1.0076 for 44 
GHz, and 1.0016 for 70 GHz. The latter is only slightly larger 
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Fig. 23. Chi2 test in pixel space, for 30,44,70 GHz. 


than the one-sigma statistical variation. For 70 GHz subsets the 
values range from 1.0017 to 1.0287. 


6.4. Pre-whitened noise 

The power of the x 2 test is limited in the sense that it collapses 
a complicated noise structure into a single number, which is not 
constrained from below. It is easy to conceive a situation where 
two errors cancel to produce a x 2 result close to one. For in¬ 
stance, a simple overestimation of the noise rms in TOI domain 
leads to an artificially lowx' 2 result. 

We complemented the validation tests by analysis of pre¬ 
whitened noise. The types of tests we performed here require 
a real-valued data set. For this reason we can perform the test in 
pixel space only. 

Consider a pre-whitened noise map, computed as 

x = M 1 m (44) 


where m is the noise map, and M is the “square root” of the co- 
variance matrix C, for which C = MM 7 . Matrix M is can be 
constructed in multiple ways. We took it to be the Cholesky de¬ 
composition of C. If matrix C accurately describes the properties 
of the noise map m, x will be a vector of uncorrelated Gaussian 
noise with unit variance, (xx T ) = I. 

We build histograms of values of x and compare them to the¬ 
oretical Gaussian distribution, in a fashion similar to the analysis 
presented in Stompor et al. ( 2002| >. We combine the 100 noise re¬ 
alizations into the same plot to improve statistics. The variance 
of the distribution yields the usual x 2 test. Comparison of the full 
distribution against the expected Gaussian distribution provides 
a more stringent test. 

We show the distributions for two selected simulation cases 
in Fig. [25] both with realistic noise. The first example is one 
where the NCVM models the noise well, the 70 GHz frequency 
map. The second example is the worst case from earlier tests, 
the FFI18/23 horn pair map with full destriping. In both cases. 
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Fig. 24. Chi2 test in pixel space, for 70 GHz subsets. Left: white 
noise. Middle : ideal and realistic noise, independent destriping. 
Right: ideal and realistic noise, full destriping. 
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Fig. 25. Distribution of noise map values before (left) and after 
(right) whitening by application of NCVM. The red curve on the 
left gives the gaussian distribution with stdev equal to that of the 
data. On the right, the red curve gives the Gaussian distribution 
with unity variance, the expected ideal result. Top: 70 GHz full 
mission, realistic simulation. Bottom: Horn pair LFI18/23, full 
destriping, realistic simulation. 


noise before whitening is very non-uniform. We show for com¬ 
parison the Gaussian distribution with variance equal to that of 
the data. The variance is dominated by the high tails of the dis¬ 
tribution, leading to a broad distribution that does not fit the 
data anywhere. On the right we show the pre-whitened noise, to¬ 
gether with the ideal Gaussian distribution with unit variance. In 
the case of 70 GHz frequency map, the distribution of whitened 
noise follows the Gaussian closely, indicating that the NCVM 
models the noise well. In the horn pair case there is a visible de¬ 
viation. There are no high tails in the distribution, the deviation 
manifests rather as a overall broadening of the distribution. In 
fact, the distribution matches very well a Gaussian shape with 
standard deviation of 1.015. 

To be more quantitative, we performed Kolmogorov- 
Smirnov test on the pre-whitened noise. We constructed the cu¬ 
mulative distribution of the values of x, and found the maxi¬ 
mum difference with respect to the Gaussian distribution. The 
probability to exceed (PTE) values computed according to 
Kolmogorov-Smirnov statistics are given in Table [6] The num¬ 
ber of degrees of freedom is N e = 921600 (100 maps of 9216 
pixels each). In general, the K-S test yields the same conclu¬ 
sions as the^ 2 test. In white noise and ideal simulation cases the 
NCVM models the noise well. The worst match is found for 70 
GHz subsets where we combine realistic simulation with full de¬ 
striping. The 18/23 horn pair again stands out as the worst case. 
The full frequency maps and submaps with independent destrip¬ 
ing are modelled well. We can identify a couple of case where 
the x 2 test indicates a better match than the K-S test. In particu¬ 
lar, for horn pair 19/22 with independent destriping and realistic 
simulation we get the excellent result x 1 - 1-0003, but only 0.01 
in the K-S test. This is likely an example of a case where two 
opposite errors compensate in the x 2 test. 
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The distribution of pre-whitened noise for horn pair 18/23, 
shown in Fig. [25] indicates that the error made in NCVM com¬ 
putation when ignoring non-idealities of realistic data, results in 
an overall underestimation of the noise, and can to some degree 
be accounted for scaling the NCVM by an appropriate constant. 
To test this hypothesis, we performed another round of K-S tests 
where we rescaled the pre-whitened noise x by its stdev value, 
bringing the variance to unity. This is equivalent to scaling the 
NCV matrix by the variance of x. The required scaling factor is 
directly given by the^ 2 result, and the^ 2 result for this modified 
test is by construction exactly one. The K-S test can still reveal 
discrepancies in the shape of the noise distribution. The PTE val¬ 
ues from the modified test are given in the last column in Table 
[6] We see that a simple rescaling of the covariance matrix leads 
to a much better agreement between the matrix and the maps. 


7. Conclusions 

We performed an analysis of noise properties of beam decon¬ 
volved maps, both in the harmonic and pixel domain. The analy¬ 
sis is built on the formalism of the artDeco deconvolution code. 
As the starting point we have a situation where a correlated 1 // 
type noise stream is first destriped to remove the correlated com¬ 
ponent as accurately as possible and then fed as input to the de¬ 
convolver. 

The main result of this work is the low-resolution noise co- 
variance matrix; this matrix describes the residual noise in the 
harmonic coefficients which is the primary output of artDeco. 
The covariance matrix models the contribution of correlated 
residual noise that is left in the data after destriping and which 
contributes strongly to the residual at lowest multipoles. The co- 
variance matrix also captures the effect of the deconvolution pro¬ 
cess itself on the noise. From the harmonic covariance matrix we 
further constructed a pixel-pixel covariance matrix, which de¬ 
scribes the noise properties of deconvolved low-resolution maps. 

As a practical application we computed covariance matri¬ 
ces for the Planck LFI four-year mission, using real LFI beams 
and noise spectra, and reconstructed detector pointings. We con¬ 
structed the matrices up to multipole f’ max =50. We constructed 
further pixel-pixel covariance matrices at resolution iV s ide=16. 

We compared the covariance products against Monte Carlo 
simulations. For validation purposes we performed a series of 
ideal simulations, which closely follow the intrinsic assumptions 
of the covariance matrix. The mean reduced x 2 values in har¬ 
monic space range from 0.9991 to 1.0071, while the one-sigma 
statistic variation is 0.0016. In pixel space the agreement is even 
better, ranging from 0.9986 to 1.0039. 

The covariance matrices are not expected to perfectly model 
realistic LFI noise due to the simplifications we have to do in ma¬ 
trix construction. We performed another series of realistic sim¬ 
ulations to find out how accurately the matrices model real LFI 
data. The mean reduced^ 2 values in harmonic space are 1.0037 
for 70 GHz, 1.0138 for 44 GHz, and 1.019 for 30 GHz. In pixel 
space the corresponding numbers are 1.0016 (70 GHz), 1.0076 
(44 GHz), and 1.01 (30 GHz). 

Various subsets of 70 GHz data gave x~ results ranging from 
1.0003 to 1.0495. We studied more closely the case with the 
largest discrepancy, and found that the main contributors to the 
discrepancy are the use of a Galactic mask, and finite baseline 
length, both related to destriping. The mask is required in LFI 
analysis to control the signal error, but our NCVM does not 
model its effect correctly. The inaccuracy arising from baseline 
length could be removed by mnning the destriping step with a 
shorter baseline. 


We considered two destriping options. In full destriping the 
full data of a frequency channel is destriped together, and a sub¬ 
set of it used for deconvolution. In independent destriping we 
use the same data subset in both processing steps. The NCVM is 
found to model the residual noise well in the latter case. The full 
destriping option yields a lower general noise level, with the side 
effect that the systematic error sources become more important. 
Consequently, the agreement between the NCVM and the data 
is poorer and the largest x 1 values are obtained with this op¬ 
tion. Fortunately, although the full destriping option yields the 
lowest residual noise for a particular submap, the independent 
destriping option is more important from the point of view of 
cosmological analysis. Noise covariance matrices are required 
input in many methods of cosmological parameter estimation. 
These methods typically input either the full frequency maps, or 
a set of submaps with independent noise, which we produced 
using the independent destriping option. 


We complemented the validation tests with Kolmogorov- 
Smirnov tests in pixel space. We used the NCV matrix to pre¬ 
whiten the simulated noise maps, and test the hypothesis that 
the whitened noise is indeed white. We observe that in cases 
where the NCVM underestimates the noise, the discrepancy is 
in most cases, at least to a first approximation, well modelled 
by an overall rescaling of the matrix. The value of the optimal 
scaling factor remains to be determined from simulations. 


It is instructive to compare the x 1 results for low-resolution 
maps to similar tests performed on the official LFI data prod¬ 
ucts, as presented in Planck Collaboration VI ( |2015[ >. The offi¬ 
cial products include pixel-pixel covariance matrices at resolu¬ 
tion V s jde=16. x 1 results from for full mission frequency maps 
are 1.042 (30 GHz), 1.027 (44 GHz), and 1.010 (70 GHz). For 
the same data sets we obtain 1.01, 1.0076, and 1.0016, respec¬ 
tively. Although our tests indicate a better agreement, the results 
must be interpreted with care. Our simulation procedure mim¬ 
ics that of official LFI pipeline. We applied same smoothing 
width (440’), the same target resolution, and the same amount 
of regularization noise. There are, however, also important dif¬ 
ferences. The results of jPlanck Collaboration VI| ( |2015| ) are 
based on 10000 realizations with a time-dependent noise model, 
while in our simulations the noise properties remains the same 
throughout the mission. A reliable comparison would require 
tests performed on the same simulation data, and NCVM con¬ 
structed with the same noise model. The FFP8 simulations used 
by Planck Collaboration VI ( 2015| ) are available in the form of 
high-resolution maps. For deconvolution, however, we need al¬ 
ternatively the TOI, or the 4D map objects constructed from it, 
which are not available 


At higher multipoles the residual noise is dominated by the 
white noise component of TOI, deformed by the deconvolution 
process. Our simulations show that deconvolution leaves rela¬ 
tively less residual correlated noise on top of the white compo¬ 
nent at high multipoles than does the ordinary binning operation. 
From this it follows that (with the exception of 30 GHz where 
the correlated residual is still large) the full noise bias at high 
multipoles can be estimated to high accuracy from the white 
noise covariance matrix, which is cheaper to compute than the 
full covariance. Also, white noise simulations are much cheaper 
computationally than full noise simulations, which require the 
destriping step. A possible future development is to equip the 
deconvolution code with an internal noise generator, which fills 
the input data structure with white noise. This would avoid the 
bottleneck of writing the simulated TOI on disk, thus allowing 
extensive MC simulations. 
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Table 6. Mean^ 2 values and PTE (probability to exceed) values from Kolmogorov-Smirnov test, comparing the pixel-space NCVM 
against 100 MC realizations. We give results for 30 GHz, 44 GHz, and 70 GHz full mission frequency maps, and for various 70 
GHz subsets (for full and independent destriping options, for pure white noise, ideal and realistic simulations). For the rightmost 
column we rescaled the NCVM by the inverse of^- 2 and re-run the K-S test for realistic simulation. The^f 2 result for this case is by 
construction exactly one. 
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30 GHz 

0.9995 
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1.0100 
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1 

0.033 

44 GHz 

0.9999 

0.602 
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0.896 

1.0076 

0.110 

1 

0.947 

70 GHz 

0.9994 
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70 GHz subsets, full destriping 
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0.271 
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1.0039 
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1 
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0.9986 
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1.0074 
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0.9986 

0.814 

0.9985 

0.914 

1.0011 

0.955 

1 

0.977 

Yr4 

1.0012 

0.108 

1.0018 

0.358 

1.0045 

0.465 

1 

0.419 

18/23 

1.0015 

0.185 

1.0022 

0.236 

1.0055 

0.076 

1 

0.368 

19/22 

0.9986 

0.065 

0.9986 

0.016 

1.0003 

0.010 

1 

0.012 

20/21 

0.9987 

0.399 

0.9991 

0.428 

1.0017 

0.180 

1 

0.332 


We computed the white noise covariance to f’ max =200. At 
higher multipoles we constructed partial matrices that cover a 
narrow multipole range around the region of interest. The noise 
bias obtained this way underestimates the noise spectrum ob¬ 
tained from Monte Carlo simulations by 1 - 3%. Should higher 
accuracy be required, the recommended approach is to estimate 
the bias through Monte Carlo simulations. That will require op¬ 
timization of the MC pipeline, as described above. 

In cases where a subset of available data is used for decon¬ 
volution, there are two options for destriping. One can either 
destripe the same data set, or the full data set. The second op¬ 
tion significantly reduces the level of residual correlated noise, 
but also leaves the different subsets correlated. We derive a co- 
variance matrix that describes the level of intercorrelation, and 
validate it on simulations involving 70 GHz horn-pair data sets. 

Finally we examined the noise properties in deconvolved 
high-resolution maps at a qualitative level. We show that the 
residual noise is correlated in the direction perpendicular to the 
beam axis. The overall noise level is low owing to the smoothing 
operation involved in map construction. 
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